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Novel behavior in electroconvection (EC) has been detected in nematic liquid crystals (NLCs) 
under the condition of comparable timescales of the director relaxation and the period of the driving 
ac voltage. The studied NLCs exhibit standard EC (s-EC) at the onset of the instability, except one 
compound in which nonstandard EC (ns-EC) has been detected. In the relevant frequency region, 
the threshold voltage for conductive s-EC bends down considerably, while for dielectric s-EC it bends 
up strongly with the decrease of the driving frequency. We show that inclusion of the flexoelectric 
effect into the theoretical description of conductive s-EC leads to quantitative agreement, while for 
dielectric s-EC a qualitative agreement is achieved. The frequency dependence of the threshold 
voltage for ns-EC strongly resembles that of the dielectric s-EC. 

PACS numbers: 47.54.-r, 61.30.Gd, 47.20.Lz 
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I. INTRODUCTION 

Electroconvection (EC) in nematic liquid crystals 
(NLCs) is an extensively studied example of electric field 
induced instabilities [ll, l2, Si- It also serves as a con- 
venient model system for investigating pattern forming 
phenomena in complex, anisotropic fluids driven out of 
equilibrium. 

NLCs have a uniaxial orientational order which is rep- 
resented by the unit vector n, the director. Consequently, 
NLCs possess direction dependent physical properties. 
The dielectric anisotropy Sa—Su— e± and the conduc- 
tivity anisotropy (Tq = (T|| — (7j_ are two key parameters in 
EC ( II denotes the value along, and _L perpendicular to 
the director). While £„ governs the electric torque acting 
on the director in the presence of electric fields, the direc- 
tor distortions lead to an elastic restoring torque exerted 
on n. Furthermore, the uniaxial symmetry of NLCs al- 
lows for a coupling between director orientation and flow, 
characterized by the 5 independent viscosity coefficients, 
resulting in a viscous torque if velocity gradients and/or 
temporal variation of n are present. 

Electroconvection has most commonly been observed 
in planarly oriented (n parallel to the substrates) layers 
of nematics with Sa < and era > 0, though it shows up 
at some other combinations of Sa, era and initial direc- 
tor orientation too [3|. In experiments usually samples 
of I0fiin< d < lOOfim thickness have been used. EC 
appears at a threshold rms value Uc of the applied ac 
voltage U of frequency / as a pattern consisting of dark 
and bright stripes with wavevector q when viewed in a 
microscope. EC patterns have great morphological rich- 
ness at the onset: normal rolls (q parallel with the initial 
n) and oblique rolls (q encloses a finite angle with n) , sta- 
tionary as well as traveling patterns have been detected 
depending on material parameters and on the driving fre- 
quency. 



Detailed theoretical studies in the past decades have 
established a firm understanding of this pattern form- 
ing phenomenon, combining the equations of nematohy- 
drodynamics (anisotropic Navier-Stokes equation for the 
flow plus the balance of torques acting on the director) 
with those of electrodynamics assuming that the NLC 
has an ohmic electrical conductivity. First it has been 
shown using a simple one-dimensional approximation 
that the EC pattern is excited by the Carr-Helfrich feed- 
back mechanism 0, Q : spatial director modulations in- 
duce a separation of space charges due to the anisotropic 
conductivity; electrostatic forces induce a flow in form 
of vortices; flow exerts a destabilizing viscous torque on 
the director against restoring elastic and electric torques. 
When this feedback becomes positive {at U > Uc) a spe- 
ciflc director fluctuation grows up; i.e., a periodic director 
modulation, characterized by a flnite tilt angle with re- 
spect to the initial alignment, develops. As the director is 
the optical axis of the NLC, the tilt leads to a modulation 
of the refractive indices. The pattern becomes optically 
detectable either with or without a polarizer due to light 
focusing/defocusing effects (shadowgraph image 0) as 
well as via birefringence if crossed polarizers are used. 
The complex physical phenomena involved in the mech- 
anism occur on three different time scales, characterized 

,2 

by the director relaxation time t^ = ^^ ^-i , the charge 

and the viscous relaxation time 



relaxation time Tq = 



o-j_ 



Tv = ^7^. Here K\\ is the splay elastic modulus, p is 
the density, 71 is the rotational, 0:4/2 the isotropic vis- 
cosity. Typically Td ^ Tq ^ Tv [i|. For example, in a 
d = 20/im sample of nematic Phase 5 (Merck & Co., Inc.) 
with a typical conductivity of crj_ = 8.2 x 10~*(r2m)~-'^ 
at T == 30°C, one has ra = 0.46s, Tq = 5.6 x lO^^s and 
Tv = 1.4 X lO-'^s. 

Though the 1-d theoretical description managed to 
capture the essence of the instability mechanism, it rep- 



resents a substantial simplification of the problem, there- 
fore the applicability of the resulting (approximate) an- 
alytical formulas are limited. This has necessitated 
later the development of a comprehensive 3-d theory, 
composed of six coupled partial differential equations 
(PDEs), known today as the standard model (SM) of 
EC [7|. The equations can be conveniently made non- 
dimensional by measuring lengths in units of cJ/tt and 
time in units of Td- Then via a linear stability analysis the 
SM can provide the frequency dependence of the thresh- 
old voltage Uc{f), that of the critical wave vector, as well 
as the spatiotemporal dependence of the director field, 
the velocity field and the charge distribution at onset. 
According to the model, a finite threshold of the primary 
instability in planar geometry requires a positive Ua and 
a negative (or slightly positive) Sa- Under these condi- 
tions, depending on the driving frequency /, two different 
solution types exist. Below the crossover frequency /c, in 
the conductive regime, the charge distribution oscillates 
with /, while the director field is stationary in leading or- 
der. In this regime, SM gives two solutions with different 
spatiotemporal symmetry (mode I and mode II), how- 
ever, for the threshold behavior only the solution with 
the lower Uc (mode I) is relevant [1,0]. Above /c, in the 
dielectric regime, the situation is reversed: the charge dis- 
tribution is stationary in leading order, while the director 
oscillates with /. Here, the model also yields two solu- 
tions (mode III and mode IV), but the threshold behavior 
is described by mode III which has lower Uc compared 
to mode IV. These electroconvection scenarios, that are 
captured by the SM, will be referred to as standard elec- 
troconvection (s-EC). Various predictions of the SM have 
been compared with experimental results obtained on a 
number of NLCs and usually an excellent agreement has 
been found - see, e.g., |3. \UX [llj . 

In s-EC, the conductive and the dielectric regimes cor- 
respond to two competing pattern forming modes which 
have different C4(/) and Uc{d) characteristics. The fre- 
quency dependence of Uc and qc has a concave (diverg- 
ing) shape in the conductive regime whereas they have 
a square root-like behavior in the dielectric one. As far 
as the thickness dependence is concerned, an inspection 
of the non-dimensional equations of the SM shows that 
by scaling only 5 of the 6 PDEs became invariant with 
respect to d; the equation for the potential still has an 
explicit d dependence via a factor Td/xq ex a±d^ - see 
e.g., equations (A6)-(A11) in [l2|. Nevertheless, in the 
conductive regime, Uc is thickness independent in the 
lowest order time Fourier approximation for the exper- 
imentally most relevant (10 — lOOfim) thickness range, 
where ra ^ Tq [7|. In contrast, in the dielectric regime 
an analogous approximation gives Uc oc d. 

At typical {d > lOfiui) cell thicknesses the conductive 
regime appears at low frequencies and the dielectric rolls 
are observable at higher /. The change of pattern type 
occurs at /c which depends significantly on a± as well as 
on d. Therefore, in order to extend the frequency range 
of the dielectric regime, thus to shift /c to lower values. 



one has to reduce either the cell thickness or the electric 
conductivity or both. When both a± and d are small 
enough, the conductive regime does not occur at all. 

The SM regards NLCs as regular dielectrics, where the 
source of the electric polarization is the electric field. In 
NLCs, however, flexoelectric polarization [l3,[lJ,[i3 may 
appear even in the absence of the electric field, if the 
director field is distorted. Flexoelectricity produces an 
additional contribution to the charge distribution and to 
the electric torque acting on the director. It is also known 
for a long time that under large enough dc voltage, flex- 
oelectricity may distort the director field in the planar 
geometry resulting in a non-convective striped pattern 
oriented parallel with the initial n [l^. The threshold 
voltage and the wavenumber of this non-convective pat- 
tern, as well as the frequency dependence of the threshold 
characteristics have been calculated [13, [l3 ■ 

The effect of flexoelectricity on EC has also been stud- 
ied and it has been shown that for the commonly used 
nematic material parameters, and for a typical sample 
thickness of d > 10/im, in the usually studied (not too 
low) frequency range of conductive s-EC, flexoelectricity 
has no significant influence on the characteristics at the 
onset of instability [I^, [20, [2l|. That gives the reason 
why the contribution of the flexoelectric polarization has 
usually been neglected in the SM. On the other hand, 
it has also been shown that when applying a dc volt- 
age, the EC threshold becomes independent of d and a± 
[l9j . Moreover, flexoelectricity leads to an appreciable 
reduction of the dc EC threshold and it also influences 
the direction of the rolls [20| . Obviously, there must be 
an ac frequency range in which the effect of flexoelectric- 
ity diminishes with the increase of the frequency. This 
frequency range, however, has not yet been studied in 
details. 

In the extended SM (i.e., with the inclusion of the 
flexoelectric polarization) two additional parameters, the 
flexoelectric coefficients ei and 63 [1, [iJl come into play. 
An analysis of the ncmatohydrodynamic equations lin- 
earized around the basic state has proven that flexoelec- 
tricity establishes a coupling between the conductive and 
dielectric modes introduced above [1, U^ resulting in a 
complex time dependence. The solutions still can be 
classified according to their spatiotemporal parity [l2| . 
The even parity solution has mode I coupled to mode 
IV (mode I-f-IV), and corresponds to the "conductive 
regime", while in the odd parity "dielectric" solution 
mode III couples to mode II (mode II+III) |22| . 

What has been said above, applies to substances with 
material parameters satisfying Eq < 0, (Xa > or Sa > 0, 
(7a < with the appropriate boundary conditions [3] . In 
contrary, in case of nematics with e^ < and ctq < 
S Ilii , where the feedback loop remains negative for all 
voltages, the Carr-Helfrich mechanism excludes the exis- 
tence of EC patterns. Nonetheless, a convection roll pat- 
tern has long ago been observed in such compounds in ac 
electric field (23l . [24| and has been reconsidered systemat- 
ically recently in a few nematics [2^, [2^, [23| . Some basic 



characteristics of these patterns, hke the orientation of 
the rolls (nearly parallel with the initial n), the contrast, 
the frequency dependence of Uc and that of the critical 
wave number have been found considerably different from 
those of the s-EC. Since this kind of pattern formation is 
not predicted by the classical SM (without the incorpo- 
ration of the flcxoelectricity), it has been called nonstan- 
dard electroconvcction (ns-EC) [23|. Recently, convec- 
tion patterns with similar properties have been detected 
in nematics built up from bent-core molecules [2^, [23, [30| 
that may exhibit giant flcxoelectricity [3 if . 

These experiments have triggered a recent reconsid- 
eration of the role of flexoelectricity within the extended 
SM for a wider range of material parameters. It has been 
shown that in some cases flexoelectricity cannot be dis- 
regarded. At high frequencies, in the dielectric regime a 
nonzero ei and 63 leads to a considerable decrease of Uc 
[3^ . What is even more important, the flexoelectric con- 
tribution to the charge separation yields flnite instability 
threshold for planar nematics with Ca < and Ea < 0, 
thus giving an explanation for ns-EC [l2| . 

Special behavior can be expected (some might go be- 
yond the applicability limit of the SM) in the parameter 
ranges where the period of the driving frequency 1// be- 
comes comparable with one of the characteristic times of 
the system Td, Tq or t^. The majority of the experimental 
and theoretical studies of EC so far apply to relatively 
thick cells of medium conductivity, thus they do not ex- 
tend to this range. We aim to fill this gap and we also 
show that flexoelectricity plays a crucial role in those 
ranges. 

In the present paper we report on experimental studies 
and numerical simulations of EC aiming to approach the 
characteristic times of the system by 1//, especially con- 
centrating on Td and in some special cases on Tq (tv is usu- 
ally very short, to approach it would require very thick 
cells and high driving frequencies, which would result in 
high threshold voltages not accessible experimentally). 
We focus on the frequency dependence of the threshold 
voltage Uc{f). The measurements have been carried out 
on five different nematic compounds: four of them ex- 
hibit s-EC, while in the fifth one, both s-EC and ns-EC 
patterns could be observed. The numerical analysis has 
been based on the extended standard model of EC which 
includes the flexoelectric effects [^, U^ ■ 

The structure of the paper is organized according to 
the pattern type in the low frequency range. After in- 
troducing the substances, the setup, and some details of 
the numerical analysis in Section [HI Section IIIII is de- 
voted to results on compounds in which conductive s-EC 
has been detected at low frequencies. In Section ITVl we 
discuss samples which exhibit dielectric s-EC at low /, 
while Section |V] describes the ns-EC instability. Finally, 
we conclude the paper with the discussion in Section IVll 



II. SUBSTANCES, EXPERIMENTAL SET-UP, 
AND DETAILS OF CALCULATIONS 



Five different NLCs have been used in the measure- 
ments. The three commercial nematic mixtures. Phase 
5, Phase 5A and Phase 4 (from Merck & Co., Inc.) 
as well as the 4-methoxy-benzylidene-4'-n-butyl-aniline 
(MBBA) have <Ta > and Eq < 0; so they exhibit s-EC 
upon application of an electric voltage. All measurements 
presented here on these compounds have been performed 
at T = 30°C. 

The fifth studied compound was the 4-n-octyloxy- 
phenyl-4-n'-heptyloxy-benzoate (8/7) [3^ which has the 
phase sequence: isotropic ~ 92 °C - nematic - 72.5 °C 
- smectic C - 62° C - crystalline. It has Eq < over 
the whole nematic temperature range, but shows a sign 
inversion from CTa < to da > as the temperature 
is increased toward the clearing point. Consequently, in 
this particular compound both ns-EC and s-EC could be 
investigated (although at different temperatures). 

Our aim is to study the EC onset behavior in the vicin- 
ity of the characteristic relaxation times of the system; in 
order to approach that we have varied the cell thickness 
and the conductivity to tune Td and Tq. 

The nematic liquid crystals have been enclosed be- 
tween two parallel glass plates coated with etched trans- 
parent indium tin oxide (ITO) electrodes. Rubbed poly- 
imide has been used to obtain a planar alignment. Cells 
in the thickness range from d ~ 3/im to d = 40/xm have 
been prepared. A spectrophotometer has been used for 
measuring d. It also enabled a systematic mapping of 
the thickness of the empty cell, thus determining its to- 
tal variation throughout the whole active area (region 
between electrodes) of the sample. Typically, we have 
found a variation of ±0.5/xm. The conductivity range 
extended from 10"9(fim)-i to 10-'^(Om)-i. The direc- 
tion of the director at the surfaces is chosen as the a;-axis. 
An electric field across the sample (along the z-axis) has 
been generated by applying a sinusoidal ac electric volt- 
age of frequency / and amplitude \^U to the electrodes. 
The cells have been placed into an oven (an Instec hot- 
stage) thcrmostatted within ±0.05°C. 

EC patterns have been studied with polarizing mi- 
croscopy using either the shadowgraph (single polarizer) 
technique or two crossed (or nearly crossed) polarizers. 
The images have been recorded with a CCD camera, dig- 
itized by a frame grabber with a resolution of at least 
768 X 576 pixels and 24 bit color depth and saved for 
further processing/analysis. 

Numerical simulations have aimed to compute the fre- 
quency dependent threshold voltage Udf) for the differ- 
ent EC patterns. The extended SM (with flexoelectricity 
included) has been used. Being interested in the onset 
behavior, linear stability analysis has been carried out as 
described in detail in [l^l . Strong anchoring of the direc- 
tor and no slip condition for the flow have been assumed 
at the bounding plates, which have been ensured by a 
Galcrkin method. The field variables have been expanded 



into sets of functions that vanish at the boundaries; the 
time periodicity has been guaranteed by Fourier expan- 
sion. Expansions in space and time have been truncated 
at the fourth and the seventh modes, respectively, which 
has been proved sufficient to obtain accuracy of better 
than 1% in C/c(/)- This has been occasionahy checked by 
increasing systcmaticahy the number of modes in space 
and time up to ten and fifteen, respectively, and moni- 
toring the changes in Uc- The growth rates ^(q, U) for 
spatial fluctuations have been calculated as a function of 
the applied voltage. The minimum of the neutral surface 
[/i(q, [/) = 0] provided the threshold voltage Uc.. The 
numerical calculations have been carried out with a For- 
tran code developed at the University of Bayreuth, using 
the same set of known material parameters (permittiv- 
ities, conductivities, clastic moduli, viscosities), as for 
earlier calculations with the SM (for Phase 5 and Phase 
5A in m, and for MBBA in 0, [H, HI]). The original 



measurements of these parameters have been reported in 
[H, M, Hi for Phase 5 and in [H, \M, EH US, SI, [H for 
MBBA. 

The splay and bend flcxoelectric coefficients, ei and 
63, should also be known, however, they are not easily 
measurable. Though various measuring techniques have 
been developed during the past decades, there is still a 
big controversy about the value and the sign of ei and 63. 
The data obtained by different methods vary consider- 
ably even for the most studied ncmatic material, MBBA: 
ei + 63 ranges from -2.3pC/m [H] to -(54 ± 10)pC/m 
4|, and 61-63 from (3.3±0.7)pC/m gl to (14±l)pC/m 
47|. Up to now, only the sign of the sum of flcxoelec- 
tric coefficients seems to be above dispute for MBBA: 
it appears to be negative (61 -1-63 < 0). For the other 
four compounds used in our experiments no data on the 
flcxoelectric coefficients are at all currently available. It 
seems though, based on measurements on a number of 
(other) nematics, that at least the order of magnitude of 
~ lOpC/m is established as a customary value for both 
61 and 63 [48| . As a consequence, during the numerical 
simulations we have regarded ei and 63 as adjustable pa- 
rameters, to be chosen to provide the best fit with the 
experiments. 



III. CONDUCTIVE S-EC AT LOW 

FREQUENCIES AND DIELECTRIC S-EC AT 

HIGH FREQUENCIES 

We start the discussion of the experimental results with 
systems exhibiting s-EC, where by choosing d and a± 
appropriately, the frequency range of the ac field covers 
1/rd and in most cases 1/Tq as well. 

The first case we show is electroconvection in Phase 
5 with d = (3.4 ± 0.5)fim and cr_L = 8.2 x 10-S(f}m)-i. 
In this case I/t^ ~ 60Hz, l/rq w 1800Hz, and 1/tv « 
2 X lO^Hz. 

Threshold measurements are presented in Fig. [TJ Be- 
low /c w 550Hz conductive s-EC (open circles) has been 



detected at onset in the form of traveling oblique rolls 
and traveling normal rolls. Above /c traveling dielectric 
s-EC (bullets) has been observed. 

A novel feature of the Uc{,f) curve has been detected 
in the low frequency range, 20Hz< / < 60Hz, where the 
threshold strongly decreases by lowering the frequency, 
in contrast to thick [d > 10/im) cells where Uc remains 
almost constant in the same frequency range. 
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FIG. 1: (Color online) Frequency dependence of the s-EC 
threshold voltage Uc measured in Phase 5 with d = (3.4 ± 
0.5)/ini - conductive s-EC (circles) and dielectric s-EC (bul- 
lets). The solid and the dashed lines are the results of the 
numerical analysis with and without the flexoelectric effect, 
respectively. The inset is the magnification of the low fre- 
quency range. The arrow indicates the crossover frequency 

/c. 

Numerical calculations have been performed with the 
Phase 5 parameter set (see Table |I]), using the flexo- 
electric coefficients 61 and 63 as fitting parameters. In 
Fig. [1] the results are plotted both without the flexo- 
electric effect (ei = 63 = 0, dashed lines, representing 
mode I and mode III) and with flexoelectric coefficients 
of 61 = — 26.5pC/m and 63 = — 23.6pC/m (solid lines, 
mode I-I-IV and mode II-I-III). For the conductive s-EC, 
the calculation without the flexoelectric effect does not 
show the bending down at low frequencies while with the 
61 and 63 values given above a good match with the ex- 
perimental data has been obtained (see the two curves in 
the inset of Fig. [T|). 

Calculations have been performed for a rf = 20/im cell 
as well, using the same parameters from Table |T] (provid- 
ing 1/rd ~ 2Hz, 1/Tq w 1800Hz, and l/r^ w 70 x lO^Hz). 
The results are plotted in Fig. [2l One can immediately 
see that the resulting Uc{f) for ei = 63 = (mode I) 
nearly coincides with the one incorporating the flexoelec- 
tric effect (mode I-I-IV) almost in the whole conductive 
frequency range, in accordance with earlier conclusions 
[19J . Nevertheless a closer look at the curves (see the 
inset of Fig. [2] with a blow up of the plot) shows that 



TABLE I: Material parameters of Phase 5 and MBBA at T 
calculations for Phase 5A. 



30° C. Phase 5 parameters have also been used in numerical 



Parameter 


Unit 


Phase 5 


MBBA 


ifii 


lO-^^N 


9.8 [34] 


5.89 [39] 


K22 


lO-i^N 


4.6 [34] 


3.7 [39] 


K33 


10"^^N 


12.7 M 


7.56 ^ 


cra/cr± 


- 


0.69 [34] 


0.5 [40] 


£± 


- 


5.25 [34] 


5.18 [40] 


Ea 


- 


-0.184 [34, 37] 


-0.48 [40] 


Qfi 


lO-'Ns/m^ 


-39 [34] 


-14.1 [43] 


Ct2 


10~^Ns/m2 


-109.3 [38] 


-80.0 [43] 


Ctz 


lO^'Ns/m" 


1.5 [38] 


-1.513 [43| 


04 


10"^Ns/m^ 


56.3 [3^ 


64.4 [43] 


as 


lO-^Ns/m" 


82.9 [38] 


57.2 [43] 


aa 


lO-^Ns/m^ 


-24.9 [38] 


-24.4 [43] 


ei 


pC/m 


-26.5 (fit) 


-14.5 (fit) 


63 


pC/m 


-23.6 (fit) 


-20.5 (fit) 



switching on the flexoelectricity results in the bending 
down of the threshold here too, similarly to the case of 
low d, only the effect occurs here at much lower frequen- 
cies. 
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FIG. 2: (Color online) Numerical calculations for the fre- 
quency dependence of the s-EC threshold voltage C/c in a 
d = 20/xm sample of Phase 5 with a^ = 8.2 x 10~*(fim)"\ 
Dashed lines are the conductive s-EC and the dielectric s-EC 
solutions with ei = 63 = 0, while the solid lines are the solu- 
tions with ei = — 26.5pC/m and 63 = — 23.6pC/m. The inset 
is the magnification of the low frequency range. The arrow 
indicates the crossover frequency /c. 

The bending down of f/c at low / when flexoelectricity 
is included can be understood from the earlier theoretical 
results [201 where it has been shown that flexoelectric- 
ity decreases the dc threshold voltage considerably (for 
MBBA the reduction is about 25%). It could also be 
proved that in this dc limit none of the dimcnsionlcss 



equations of SM depend on a \ (P , hence the dc thresh- 
old is independent of d and cr j_ [1^ . In case of ac driving 
there is a wide frequency range (in the middle of the con- 
ductive regime) where flexoelectricity has much less influ- 
ence on Uc] it leaves the ac threshold almost unaltered, at 
least for not too thin cells. This means that at increasing 
/ the contribution from the flcxoelectric charge separa- 
tion (important at dc) fades away compared to that of 
the Coulomb charge density in a narrow frequency range 
at low /. The character of this transition from dc to ac 
driving has, however, not been studied rigorously either 
experimentally or theoretically so far. 

In order to understand better the behavior of Uc we 
have carried out further simulations, carefully analyzing 
the low frequency range for a series of thicknesses when 
flexoelectricity is included. Results are shown in Fig. [3] 
for three thicknesses. The figure indicates that the tran- 
sition from ac to dc (the bending down of Uc to its dc 
value) occurs for any thickness as expected; only the ef- 
fect is compressed into a very narrow frequency range for 
large thicknesses. The corresponding values of \/t^ are 
also indicated in Fig. [3] (circles). Their values fall into 
those frequency ranges (different for each d) in which the 
strong bending down of Uc occurs. We note, that accord- 
ing to its definition, rd is the relaxation time of a homo- 
geneously deformed director state (wavenumber |q| =0). 
The relaxation time of the pattern, l/rp, is shorter, be- 
cause the growth/decay rate (/j, = I/tjj) of a periodic 
(|q| ^ 0) distortion increases with |q| [43, [50|. Therefore 
1/Tp > 1/Td would be a more appropriate characteristic 
frequency of the system than l/rd. A rigorous quantita- 
tive comparison of l/rp with / is, however, problematic 
since |q| and hence Tp are frequency dependent. 

The crossover between the regimes where the Coulomb 
charge separation or the flcxoelectric charge separation 
dominates, is reflected in the curvature of the C4(/) 
curves. It is convenient to use the location of the in- 
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FIG. 3: (Color online) Numerically calculated frequency de- 
pendence of the s-EC threshold voltage C/c for samples of var- 



with 



;.2 X iQ-^'inmy 



ious thicknesses in Phase 5 
ei — — 26.5pC/m and ea = — 23.6pC/m. The circles and 
the star symbols on the curves denote 1/rd and the inflection 
point at /infl, respectively. Inset: the thickness dependence of 

l//infl. 



flection point at find (where the curvature changes sign, 
indicated by stars in Fig. [3]) as a quantitative measure. 
As one sees values of /infl are close to l/rd, but interest- 
ingly, l//infl is linear with d (see in the inset of Fig. [3]), 
in contrast to the (P dependence of Td. 

Finally, we mention that for the high frequency end 
of the conductive regime the influence of flexoelectric- 
ity becomes again noticeable though less dramatic. It 
considerably enhances the d-dependence (see in Fig. [3]) 
which would be almost negligible at ei = 63 = 0, espe- 
cially for higher thicknesses. Furthermore it shifts /c to 
lower frequencies as seen in Fig. [TJ 

The shift of /c is mainly due to the fact that in 
the dielectric regime of s-EC the flexoelectric effect de- 
creases the threshold voltage drastically - compare the 
dashed curve (ei = 63 = 0, mode III) and the solid line 
(ei = -26.5pC/m and 63 = -23.6pC/m, mode Il-flll) 
for the dielectric solutions above /c in Figs. [T] and [2] The 
reduction reaches about 16% even in the relatively thick 
(d = 20/xm ) sa mple, which agrees with recent findings 
for MBBA [34|. Consequently, flexoelectric effects can- 
not be neglected in case of dielectric s-EC even in the 
thick samples. We also mention here that no peculiar ef- 
fect on Uc{f) is observed in the vicinity of 1/rq w ISOOHz 
(see Fig. III). 

Summarizing the results above one can conclude that 
with the above mentioned values of a±, ei and 63, an 
excellent agreement has been found between the experi- 
mental results (symbols) and the numerical calculations 
(solid lines in Fig. [1]) for Phase 5 in the whole experi- 



mentally covered frequency range. Note that the fitted 
ei and 63 values are realistic in a sense that they are of 
the same order of magnitude as those measured for other 
nematics. 

Similar results have also been obtained on Phase 5A 
and MBBA samples in which only the conductive regime 
of s-EC has been detected. Phase 5A has been investi- 
gated in a sample of d = 3.1/im. At dc voltage as well as 
at very low frequencies (up to / w SHz) static stripe pat- 
terns, the so called flexoelectric domains [13, ^M have 
been observed at a threshold U[. Above / = 3Hz 
only conductive s-EC patterns have been seen: up to 
/ — 40Hz stationary oblique rolls and above / ~ 40Hz 
traveling oblique rolls and traveling normal rolls. The 
measured frequency dependence of the threshold voltage 
is presented in Fig. [Hby open circles. Below / « 50Hz 
it appears that Uc decreases considerably as approaching 
/ = (see the inset in Fig. S]), in a similar fashion as 
seen for the Phase 5 in Fig. [TJ Again the bending down 
occurs in the frequency range of l/rj. 
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FIG. 4: (Color online) Frequency dependence of the threshold 
voltage C/c in the conductive regime for a sample of Phase 5A 
with d — S.lfim. Circles stand for experimental data. The 
solid line represents the Uc{f) curve calculated numerically 
with Phase 5 parameters (same as in Figure 1, with nonzero 
ei and 63) and with a± = 9.2 x 10~*(f2m)~^. The dashed 
line is the calculated solution for ei = £3 = 0. The inset is 
the magnification of the low frequency range. 

Phase 5A has the same chemical composition as Phase 
5 except that Phase 5A is doped with some ionic salt 
in order to increase the electric conductivity. Therefore, 
in the numerical calculations for Phase 5A we have used 
the material parameters of Phase 5 (given in Table |I]) 
with a± = 9.2 x 10^^(Om)~^, and the flexoelectric coef- 
ficients ei = — 26.5pC/m and 63 = — 23.6pC/m that have 
been obtained from the best fit in Fig. [TJ The calculated 
threshold curve for conductive s-EC (solid line in Fig. ^ 
agrees very well with the experimental data in the whole 
frequency range. 



For reference, in Fig.|4]we have also plotted the conduc- 
tive s-EC threshold Udf) without the flexoelectric effect 
(i.e.; for ei = 63 ~ 0; dashed line). The same conclusions 
can be drawn here as those for Phase 5: the Udf) curve 
with ei = 63 = deviates significantly and systemati- 
cally from that with nonzero ei and 63 (and also from 
the experimental data), especially in the low frequency 
range (below / w 50Hz) and at high frequencies (above 
600Hz). 

Up to now the most studied nematic EC material is 
presumably MBBA which we have investigated too, in a 
cell of d = 3.2/im. At dc voltage neither flexoelectric do- 
mains, nor periodic EC patterns have been detected up to 
U = 80V. However, at and above / = O.lHz a conductive 
s-EC pattern has emerged. Frequency dependence of the 
threshold voltage measured in this thin layer of MBBA is 
shown in Fig. [5] (circles) . Again, below / w 30Hz Uc{f) 
bends down considerably while approaching / ^ 0, in a 
similar fashion as seen for Phase 5 and Phase 5A. 
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FIG. 5: Frequency dependence of the conductive s-EC thresh- 
old voltage Uc measured on a sample of MBBA with d = 
3.2/^ni (circles). The solid line represents the Uc{f) curve 
calculated numerically with MBBA parameters from Table 
U and with a± = 9 x 10"**(nm)"\ ei = -14.5pC/m and 
63 — — 20.5pC/m. The dashed line is the solution for the 
same parameter set, but with ei = es = 0. 



The MBBA parameter set from Table H] has been used 
in the numerical calculations. The solid line in Fig. [5] 
obtained by the simulations for ei = — 14.5pC/m, 63 ~ 
— 20.5pC/m and a± = 9x 10~®(fim)^^ represents a fairly 
good match with the experimental data, in contrast to 
the curve with ei = 63 = (dashed line in Fig. [5]) which 
deviates significantly. One has to note here that the fitted 
values of ei and 63 fall in the mid-range of measured 
values for MBBA [i^ , and in particular agree reasonably 
well with the values |ei| ~ 15pC/m and |e3| ~ 30pC/m 
reported in [53, \EM ■ 



IV. DIELECTRIC S-EC AT LOW 
FREQUENCIES 

With a careful preparation procedure we were able to 
produce thin (d = 3.4/im) Phase 5 samples with a con- 
ductivity as low as about <t± Ki 1 x 10~*(r2m)~^. These 
samples exhibited convection-free fiexoelcctric domains 
at a threshold U{ at dc and at very low frequency (up 
to / = 7.5Hz) ac driving, similarly to the situation in 
Phase 5A. Above / = 7.5Hz, in a narrow frequency range 
(up to / « 20Hz) fiexoelcctric domains and traveling di- 
electric s-EC oblique rolls coexist with slightly different 
thresholds. At higher frequencies dielectric s-EC could 
be detected in the form of traveling oblique, or traveling 
normal rolls. 

The frequency dependence of Uc in this sample (de- 
picted with bullets in Fig.lS]) is, however, rather unusual 
and surprising, considering that both the SM and pre- 
vious experiments on dielectric s-EC provided a square 
root-like Udf) function. Instead, here Uc{f) bends up 
as / ^ 0, i.e. it is a non-monotonic function with an ex- 
pressed minimum at /min in the frequency range of l/rj. 
Moreover, at higher frequencies (above the minimum) 
Uc{f) is linear within the experimental error. Perform- 
ing numerical simulations with the Phase 5 parameter 
set in Table H] and with (t± = 1 x 10^®(r2m)~^, we have 
found that the calculated dielectric threshold curve Udf) 
(shown as the solid line with the lowest Uc in Fig. [5]) is 
monotonic, has a smaller slope (above 60Hz) and pro- 
vides lower thresholds (especially at high /) than found 
experimentally. 
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FIG. 6: (Color online) Frequency dependence of the dielec- 
tric threshold voltage Uc measured in a d = 3.4/xm sample of 
Phase 5 (bullets). Lines represent results of numerical calcu- 
lations for Phase 5 with ei — — 26.5pC/m, 63 = — 23.6pC/m, 
and a± = 1 X 10~*(nm)~^ for a series of aa/o-± values. 



In order to test whether the extended SM is able at all 
to provide an Uc{f) dependence with at least a qualita- 
tive similarity to the measured data, we have performed 



simulation experiments. First we have checked the influ- 
ence of the conductivity systematically in a wide range 
of 10"^2(f^nj)-i <a^< 10~S(17m)-\ where only dielec- 
tric s-EC is expected. We could verify that the calculated 
Ucif) remains monotonic in this a± range, and dUc/df 
is not affected significantly. Consequently, the novel fre- 
quency dependence of the threshold (bending up) cannot 
be interpreted by a pure reduction of the conductivity. 

We have also varied the anisotropy of the conductiv- 
ity motivated by recent numerical results |12l | indicating 
that a decrease of (Ta/cj. increases the slope dUc/df and 
tends to make Uc{f) linear. In Fig. [6] we present the cal- 
culated threshold curves for five additional (Ta/a± values, 
leaving all other parameters in Table |l] unaltered. It is 
seen that Uc{f) as well as dUc/df shift to higher values 
with the decrease of <7a/<^j_, moreover, Uc{f) becomes 
non-monotonic for aa/o'± < 0.2 with a bending up of 
Uc{f) at low frequencies. Note that the curve calculated 
for aa/cF± = 0.1 reproduces qualitatively all features of 
the experimental Uc{f), except that it is shifted to higher 
voltages. 

To demonstrate the role of flexoelectricity in the novel 
frequency dependence of the dielectric regime, in Fig. [7] 
we have plotted Udf) calculated with (Ja/o'± = 0.1 both 
for mode III (i.e., without fiexoelectricity, dashed line) 
and for mode II-I-III (i.e., with flexoelectricity, sofid line). 
Obviously, when flexoelectric effects are neglected, the 
threshold is even higher, Uc{f) decreases monotonously 
when / — ^ and the bending up is not reproduced, so 
there is no agreement with the experimental data even on 
a qualitative level. On the other hand, if one takes into 
account the flexoelectric coefficients given in Table U an 
acceptable qualitative agreement is achieved: the curve 
has a minimum at about the same /min frequency where 
the experimental data, and for / > /min the steepness 
dUc/df is roughly the same as that in the experiments. 
Therefore, the bending up can undoubtedly be attributed 
to flexoelectricity, and /min is a natural indicator for the 
frequency range where the crossover between dominantly 
flexoelectric and dominantly Coulomb charge separation 
mechanisms occur. 

Finally we mention that (similarly to the case discussed 
for Fig. [5]) no peculiar effect of flexoelectricity on Uc{f) 
is observed in the vicinity of l/rq « 220Hz (see Figs. [6] 
andEl). 

Systematic numerical calculations at (7a/ a± = 0.1 have 
shown that while keeping aj_ (and thus Tq) constant, /min 
shifts to higher / with the decrease of d (i.e., of Td) - 
see Fig. [H Despite of the fact that l//min increases 
linearly with d up to d = 8/im (similarly to l//infl of 
conductive s-EC), one finds the values of 1/ fmin{d) close 
to Td(d), and the best fit l//min oc d^ does not differ 
significantly from Td{d) (see inset of Fig. [8]). At d — 20/im 
Uc{f) becomes a monotonically increasing function (at 
least above / = SHz). One has to mention here that 
in our calculations we have not resolved the very low 
frequency range (^^ IHz), where the theoretical analysis 
becomes more difficult due to numerical problems and the 
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FIG. 7: (Color online) Numerical calculations of the threshold 
Ucif) for Phase 5 with aa/a_L = 0.1, ax = 1 x 10"*(nm)-\ 
and d — 3.4/im: dashed line is dielectric s-EC with ei = ea = 
0; solid line is the dielectric mode with ei = — 26.5pC/m, 
63 = -23.6pC/m. 



experiments also require special care (and are, therefore, 
typically avoided) . From Fig. [8] one also sees that the 
steepness of Udf) in the high frequency range (where 
the function is close to linear) also depends on d; dUc/df 
increases with d. 
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FIG. 8: (Color online) Frequency dependence of the dielec- 
tric s-EC threshold, calculated with Phase 5 parameters for 
o'a/o'± = 0.1 and with ei = — 26.5pC/m, 63 = — 23.6pC/m, 
CTx = 1 X 10~*(Slm)~^ for different sample thicknesses d. Cir- 
cles denote 1/rd for each d. In the inset the thickness depen- 
dence of l//min (bullets) and that of rd (solid line) is presented 
together with the best fit l//min oc d^ (dashed line). 

The other situation is when Tq is tuned via a± while 
keeping d and thus Td constant (l/rd = 75.5Hz). The 
infiuence of (J± on Udf) is illustrated in Fig. [5] With 



the increase of the conductivity (i.e, of 1/Tq), both /„ 
and the lowest threshold value C/c(/min) increase. 
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FIG. 9: (Color online) Results of numerical calculations 
for Phase 5 with aa/o-± — 0.1, ei = — 26.5pC/ni, 63 — 
— 23.6pC/ni and d = 3.4^ni for three different conductivities 
a-±. 

Samples of the nematic Phase 4 have also been inves- 
tigated at r = 30° C. Measurements on the frequency 
dependence of Uc are shown in Fig. [TUlfor a rf = 3.2/im 
cell (circles). At dc voltage (solid triangle) and at low 
frequencies (up to f — 15Hz) flexoelectric domains [see 
in Fig. [TTT a)] have been observed with a tendency of Uf 
increasing with the frequency. At / = 15Hz the thresh- 
olds for flexoelectric domains and for traveling dielectric 
oblique rolls almost coincide (see the coexisting patterns 
in Fig. [TlT b) - a similar behavior has also been reported 
in [5J|). At higher / flexoelectric domains could not be 
detected. Instead, traveling dielectric oblique rolls [Fig. 
Illf c)] or traveling dielectric normal rolls have been ob- 
served. The threshold curve for these traveling rolls has 
a minimum similarly to that observed in dielectric s-EC 
of Phase 5. In contrast to the thin cell, however, ex- 
perimental data obtained on a d = 12.4/xm thick sample 
of Phase 4 (star symbols in Fig. [TU] and in its inset) 
demonstrate a monotonically increasing Uc{f) (at least 
for / > lOHz). At the same time, Fig. [TOl represents 
an experimental support of numerical results shown in 
Fig. [S] obtained though with Phase 5 parameters (unfor- 
tunately, in Phase 5 samples with d > lOfim, dielectric 
s-EC could not be obtained at low enough /). 



V. NONSTANDARD EC 

Previous measurements have shown that in samples 
of 8/7 with typical thickness (d > lOfim), ns-EC takes 
place in the whole investigated frequency range at T < 
85°C, characterized by a linear Uc{f) dependence. In the 
temperature range 85°C< T < 90° C both ns-EC and s- 
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FIG. 10: (Color online) Frequency dependence of the s-EC 
threshold voltage Uc measured in a sample of Phase 4 with 
d = 3.2/im (circles). At dc voltage, the threshold for the 
flexoelectric domains is also indicated (solid triangle). For 
comparison, experimental data obtained on another sample 
of Phase 4 with d — 12.4/im is also plotted (stars). The inset 
displays experimental results on the d = 12.4/im cell in the 
whole frequency range studied. 




FIG. 11: Snapshots of patterns taken slightly above the 
threshold in a d = 3.2/^m sample of Phase 4. (a.) flexo- 
electric domains at dc voltage; (b.) coexisting flexoelectric 
domains and traveling dielectric oblique rolls at / = 15Hz; 
(c.) traveling dielectric oblique rolls at / = 70Hz. The arrow 
denotes the initial director. The physical size of the images is 
(65 X 65)iim^. 



EC could be observed depending on / and d, while above 
T = 90°C always s-EC has appeared at onset |26| . 

In accordance with the above experimental findings, in 
our d — 3.4/^m sample of 8/7 ns-EC has been found be- 
low T w 90°C. At higher frequencies, similarly to thick 
cells, a linear frequency dependence of the threshold volt- 
age has been detected which is illustrated in Fig. [T^] for 
three different temperatures. However, at lower / the 
frequency dependence of Uc is different from the one ob- 
served in thicker samples [2^. The threshold voltage 
increases abruptly here when / — > as shown in Fig. 
1121 i.e., Ucif) becomes a non-monotonic function with a 
minimum fJcmin at /min, similarly to the case of dielectric 
s-EC in Phase 5 and Phase 4 samples. 

Studies on the temperature dependence of Udf) in the 
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FIG. 12: (Color online) Frequency dependence of the ns-EC 
threshold voltage Uc measured in a sample of 8/7 with d = 
3.4/^m for three different temperatures T. Solid lines represent 
linear fit to data at high frequencies (above the minima of the 
curves) . 



sample of 8/7 with d = 3.4yum have revealed that: (i.) 
f^cmin reduces with the increase of T; (ii.) /min increases 
with T; (iii.) the steepness of Uc{f) in the high frequency 
range (where the function is linear) is also temperature 
dependent; dUc/df diminishes linearly with the increase 

of r. 



nomcnon of the ns-EC pattern at Uc in thicker samples, 
in Fig. [13] we compare this C/js with Uc of ns-EC mea- 
sured in samples of different thicknesses at T = 76° C. 
It is seen that the novel non-monotonic behavior is re- 
stricted to the two smallest thicknesses (d = 2.4/im and 
d — 3.4/im). /mill increases as d is lowered (it is about 
220Hz for the d ~ 3.4/im sample, and is roughly 500Hz 
for d = 2.4/im). Therefore, one can conclude that /min 
depends stronger on d, than on the material parameters 
(compare with the much weaker temperature dependence 
of this minimum in Fig. I12p . It can be noted that while 
C^cmin and /min grow, dUc/df decreases with lowering the 
sample thickness d. 

The temperature range where ns-EC is observable ex- 
tends to higher T in thin samples compared to thicker 
ones. To demonstrate this, we plotted in Fig. [H] the fre- 
quency dependence of the threshold voltage at T w 90° C 
for three different sample thicknesses. In the thick, d = 
40/im, cell s-EC exists in the whole accessible frequency 
range. At an intermediate thickness of 13.2/^m a s-EC 
to ns-EC transition has been observed at / w 500Hz. 
In the thin, d = 3.4^m, cell ns-EC has been detected 
in the whole frequency range. In this latter case a non- 
monotonic Uc{f) is observed, just as at lower tempera- 
tures, or in the dielectric s-EC of Phase 5 and Phase 4. 
Fig. [HI demonstrates that a s-EC to ns-EC transition can 
also be induced by reducing the sample thickness (while 
keeping the temperature and frequency fixed), in addi- 
tion to the already reported temperature, frequency and 
electric field induced transitions |26j]. 
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FIG. 13: (Color online) ns-EC threshold curves Uc{f) mea- 
sured in samples of 8/7 of different thicknesses d at T = 76°C. 
Solid lines represent linear fit to data. 

In an even thinner {d = 2.4/im) sample of 8/7 no peri- 
odic EC structure could be observed, instead, a large- 
scale, non-periodic pattern has appeared at a critical 
voltage Uis- Since a similar non-periodic, large-scale 
pattern has also been observed as a concomitant phe- 



FIG. 14: Frequency dependence of the threshold voltage Uc 
measured in 8/7 samples of different thickness {d = 3.4/im, 
13.2/im and 40/im) at T f» 90°C (in the vicinity of the s-EC to 
ns-EC transition). Open symbols denote the s-EC threshold, 
while closed symbols stand for ns-EC. 

Recent theoretical advances (the analysis of the ex- 
tended SM) have shown that the flexoelectric charge 
separation mechanism could be responsible for the oc- 
currence of ns-EC [131 . Numerical simulations based on 
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the extended SM could actually yield a good agreement 
with experimental results on a thick sample of the com- 
pound 8/7 using a reasonable guess for unknown ma- 
terial parameters. Extension of the numerical analysis 
to smaller thicknesses and lower frequencies is, however, 
not so straightforward in the parameter range of ns-EC as 
for s-EC. In addition some important parameters (e.g., 
viscosities and flexoclcctric coefficients) should first be 
determined by independent measurements. As a conse- 
quence a theoretical f/c(/) dependence for Fig. [Til could 
not yet been provided; it remains a target for further 
studies. 



VI. DISCUSSION 

We have reported about measurements of the threshold 
voltage of electroconvection patterns in various nematic 
liquid crystals focusing mainly on thin (d « 3.1 — 3.5/im) 
samples in which I/t^ is relatively large. These samples 
have exhibited patterns of various types in the vicinity 
of l/Ttj: conductive s-EC, dielectric s-EC as well as non- 
standard EC. For all compounds and pattern types we 
have found a novel frequency dependence of the threshold 
voltage at low frequencies, which deviates significantly 
from that typically seen in thicker cells with smaller I/t^. 
This indicates that this novel phenomenon is quite ro- 
bust. The character of the novel Udf) behavior seems, 
however, to depend on the type of the pattern. On the 
one hand, whenever conductive s-EC developes at onset, 
Uc{f) expresses a pronounced decrease (bending down) 
if the frequency is lowered, thus Uc{f) has an inflection 
point. On the other hand, if the onset pattern is dielec- 
tric s-EC in the vicinity of l/rd, Uc{f) exhibits an abrupt 
increase (bending up) towards lower /, so Uc{f) has a 
minimum. In case of ns-EC, the low frequency behav- 
ior is similar to that of the dielectric s-EC. As according 
to recent theoretical results the ns-EC is actually a di- 
electric mode [12[, one can conclude that whether Udf) 
bends down or up at / -^- 0, depends on the spatiotem- 
poral symmetry of the solution at onset. The significant 
changes in the curvature of Udf) occur in the frequency 
range of l/rd, even though some characteristic frequen- 
cies (/infl for the conductive or /min for the dielectric 
pattern) do not scale precisely with t^- 

Besides the experiments, the threshold voltages have 
also been calculated numerically for s-EC in a linear sta- 
bility analysis using both the standard model of electro- 
convection (neglecting flexoelectricity) , and the extended 
SM (including the flexoelectric effect). The simulations 
have shown that in the absence of flexoelectricity thin 
as well as thick cells exhibit the "regular" Udf) behav- 
ior, without the novel dependence at low /. Including, 
however, flexoelectricity with properly adjusted flexoelec- 
tric coefficients the pronounced low frequency decrease of 
the conductive s-EC threshold could be reproduced in a 
quantitative agreement with the experimental data. 

Unfortunately, a similar quantitative match could not 



be achieved for dielectric s-EC patterns. Nevertheless, 
simulation experiments have shown that even the abrupt 
increase of the dielectric Udf) in thin cells at low / can 
be qualitatively reproduced by the extended SM if one 
uses (together with the same flexoelectric coefficients as 
for the conductive s-EC) a (Ta/crj^ value which is lower 
than expected from the measurements. One has to em- 
phasize here that the fitted ei and 63 parameters fall into 
the range of values measured for nematics, however, we 
do not claim that the fitted values are unique, some other 
ei, 63 combinations might provide almost as good match 
with experimental data. Furthermore, in the numerical 
calculations only a subset of parameters (ffa, cr_L, ei, 63) 
have been varied. One can not exclude that by slight 
adjustment of the other parameters (allowed by experi- 
mental errors in measuring the values in Table |T| a better 
agreement between measurements and calculations could 
be obtained. 

The numerical simulations presented above have con- 
firmed that the experimentally detected novel low fre- 
quency dependencies of Udf) a-re consequences of the 
flexoelectricity and represent a competition between flex- 
oelectric and Coulomb charge separation mechanisms. 
The frequency range in which it occurs is in the order 
of 1/Td, thus it can be tuned by the proper choice of the 
sample thickness d. 

In the paper we have focused on the EC threshold volt- 
ages only. For a complete analysis of the onset behavior 
one should also confront the experimental and theoretical 
data regarding the critical wave- vector qc of the pattern. 
The special experimental conditions, namely the small 
wavelength of the pattern (about the same - few microns 
- for conductive and dielectric s-EC as well as for ns- 
EC), the fast traveling speed of the rolls, the relatively 
weak contrast (weaker than in thick samples) unfortu- 
nately did not allow to obtain qc values with required 
precision by polarizing microscopy. 

From the misfit of the dielectric thresholds one can 
draw a double conclusion. On the one hand, the obtained 
qualitative agreement may indicate that the extended SM 
does contain the main ingredients of the physical mecha- 
nism responsible for the formation of dielectric patterns 
even under our experimental conditions. On the other 
hand, however, the lack of quantitative match implies 
that the applicability limits of the model could be ap- 
proached by testing thin cells at low frequencies and in- 
clusion of additional phenomena (resulting in an appar- 
ent reduction of Ca/c^) would be needed for a complete 
description. One such candidate could be the (unipolar) 
charge injection at the electrodes which is not included 
in the SM. For a low value of the overall conductivity (as 
in our case for dielectric s-EC) the weak charge focusing 
effect may be disrupted by charge injection at dc or at 
low frequency ac fields as it has been suggested recently 

We also have to note here that in our thin samples for 
all compounds traveling EC rolls have been detected at 
the onset. It has been proven that traveling waves can 
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be explained if the assumption of Ohmic conductivity is 
given up and the SM is replaced by the weak electrolyte 
model (WEM) [5a| taking into account ionic phenom- 
ena. It has also been shown that for typical EC samples 
{d > 10/xm) in the conductive s-EC regime the Uc and Qc 
values provided by the WEM differ by less than 1% from 
those calculated with the SM 
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However, in this com- 
parison flexoclcctricity has not been included. Neither 
has the test been performed for smaller d values, nor for 
dielectric s-EC. Therefore, a more complete description 
of the observed phenomena would need the inclusion of 
the flexoelectric effect into the WEM, and an extension 
of WEM based numerical studies to dielectric s-EC as 
well as to ns-EC, which is going far beyond the scope of 
the present paper. Such simulations would require some 
knowledge about additional material parameters, such as 
ionic mobility, recombination rate, equilibrium ion den- 
sity. In this respect the measurements of the dc and low 



frequency ac conductivities of thin nematic layers could 
be important to clarify the impact of weak electrolyte 
effects. Future studies in these directions could possibly 
bridge the quantitative gap between the experiments and 
numerical data presented in this paper for the dielectric 
s-EC. 
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